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Abstract 

The tunneling splitting of the energy levels of a ferromagnetic particle in the 
presence of an applied magnetic field - previously derived only for the ground state 
with the path integral method - is obtained in a simple way from Schrodinger theory. 
The origin of the factors entering the result is clearly understood, in particular the 
effect of the asymmetry of the barriers of the potential. The method should appeal 
particularly to experimentalists searching for evidence of macroscopic spin tunneling. 
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In spite of early speculation ]I[ that the one loop approximation of a path integral 
calculation in quantum mechanics would not yield exactly the same result as a WKB 
semiclassical calculation, it is now generally agreed - and has indeed been verified 
by application to important examples - that the two approaches yield exactly the 
same result in the approximation of correct linear matching of WKB solutions across 
turning points. Thus with this proviso the latter method was shown || to lead 
to the same level shifts in the cases of the double-well and periodic potentials as 
the path integral method with expansion around an instanton. The effect of the 
nonlinearity which gives rise to instanton or other classical solitonic configurations 
arises in Schrodinger quantum mechanics through the boundary conditions imposed 
on the perturbatively calculated wave functions. It is these boundary conditions 
therefore which generate the nonperturbative effects characteristic also of tunneling. 

Once the Hamiltonian of a theory with tunneling effects has been established, 
it is still a long way to obtain the level shift or decay rate with the path integral 
method. Such calculations are very instructive, particularly as models for more com- 
plicated field theoretical contexts, but can lack transparency in other contexts where 
a confrontation with experiments is anticipated. Thus recently long and complicated 
instanton path integral calculations || ||, ||, [| |7], ||, §, |H| |TT| were performed in 
order to derive the tunneling splitting of the ground state of a magnetic particle in 
a crystal field with or without an applied magnetic field. This is a case which is of 
considerable experimental interest since the tunneling rate should be directly exper- 



imentally measurable in resonance experiments[ll]. It seems it is not so well known 



that the calculation of such tunneling effects can also be performed in the context of 
Schrodinger theory with appropriately chosen (in this case periodic) boundary con- 
ditions, and in fact much more directly so by comparison with established cases as 
demonstratred recently in the case with no applied magnetic field[0. On the basis 
of our experience with both types of methods (for path integral calculations see also 
refs. [ 13, |T5|), we expect that any calculation which is manageable in the path 
integral method can also be dealt with in Schrodinger theory and vice versa. In fact, 
the latter even has certain advantages in yielding with equivalent ease (or complex- 
ity) the corresponding results for excited states whereas in the path integral method 



this would necessitate the consideration of periodic or nonvacuum instantons[|14], [15 
which implies a further degree of considerable calculational complication. 

Quantum spin tunneling (or more generally tunneling for internal degrees of free- 
dom) has only recently been realised to be a fascinating quantum mechanical phe- 
nomenon. In the following we demonstrate that the results obtained in refs.[]3], |J 
and [0] for spin tunneling between classically degenerate states in the presence of 
an applied magnetic field can essentially be obtained in a very simple way by com- 
parison of the Schrodinger equation with that of a periodic potential for which the 



level splitting is well known, i.e. the Mathieu equation|16|. Apart from supplying a 



check of the path integral result, this type of derivation should appeal particularly to 
experimentalists. Moreover it has the enormous advantage of yielding automatically 
the extension to excited states which - as remarked earlier - is much more difficult 
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to achieve with the path integral method. The calculation also gives a clearer picture 
of which parameters have to be large and which small for the results to have validity 
in the sense of asymptotic expansions. 

In the following we concentrate on the case of a magnetic particle (e.g. an ion) 
embedded in the field of a crystal and a weak applied magnetic field h = qhb-B 
perpendicular to the easy axis X. As explained in refs. || £| for a specific model 
(among others) the appropriate Hamiltonian is 

H = -AS 2 X + BS 2 Z - hS y (1) 

where A and B are positive anisotropy constants. One should note that two anisotropy 
axes are needed to produce tunneling when ft = 0. After transforming the spin op- 
erators Si to canonical operators p, <f) [|17[] the relevant semiclassically approximated 
Hamiltonian can be written 
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where s is the spin of the particle which is assumed to be not too small. The classical 
phase space is the unit sphere. The potential (which is two-fold degenerate for ft = 
in the domain < < tc with easy axis X, i.e. the degenerate minima at and n) 
can be seen to result from the Hamiltonian H obtained by replacing Si by sei(9, 4>), 
ej a unit vector expressed in polar coordinates, and replacing 9 by its classical ground 
state value |. The derivation of the field dependence of the mass (which does not 
concern us here) is described in refs.|| |J. In the following we set 

B ft 
o = — , a = ; (5) 
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and consider the experimentally interesting case of b > 1. Thus the mass m(<p) is to a 
first approximation given by We demonstrated earlier |12| that corrections to this 
are small and can be taken into account in terms of elliptic functions (a modification 
which could, in principle, also be considered here but would overload this note with 
algebraic technicalities). We ignore these corrections in the first place. 

The quantum states of the system (i.e. the degenerate states separated by in- 
finitely high barriers) are determined by the oscillator approximation of the system 

around minimum positions 0o of V given by (with s(s + 1) ~ (^s + ~ 
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and so 



y 0,2n 



2nn + arcsin a, 4>o,2n+i = (2n + l)7r — arcsin a, n — 0, 1, 2, 



(7) 



with 
at which 



COS0 O ,2n = +Vl - a 2 , COS 00. 



2n+l 



V"(0 O ) = 2As(s + 1) [1 - a 2 ] = 8h 2 m .B (8) 

In the problem at hand < < 2tt and there are only two minima as in the case of 
the well known double well potential. The Schrodinger equation = defined 
by H is then approximately 

~E - V(<f)o 
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* ~ 



This determines immediately the oscillator approximated eigenvalues as 

h 2 



E. 



(0) 

2n+l 
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(9) 



(10) 



in agreement with the "semiclassical ground state energy" (n = 0) given in ref.|| 
(last formula) except that our B (assumed to be much larger than A) is there B + A. 
This perturbation theoretical expression ignores tunneling. 

A typical aspect of any tunneling formula is the exponential of the euclidean 
action of the classical vacuum pseudoparticle (which tunnels through the barrier). 
This factor supplies the classical approximation of the transition amplitude equivalent 
to the wave function approximation given by the WKB exponential in quantum 
mechanics and must be such that it vanishes (i.e. the level splitting) in the limit of 
infinitely high barriers. In the present case this implies that l? m has to be large and 
thus s or As(s + 1). The argument of the exponential must therefore contain s. In 
the path integral results of refs.|3], [| the factors raised to powers of 2 (s + ~) play 
this role. In fact the factor in eq. (9a) of ref.J| can be shown to approximate the 

WKB factor exp(— 2 ^^ 2 ^ ) in the limit of vanishing magnetic field h (one may note 
that a = 1 determines a critical value of h beyond which the degeneracy of minima 
for h ^ — cf . eq.(D) - is removed as discussed by Schilling |11]). These observations 

(i.e. the lifting of the degeneracy of oscillator levels) 



cf. eq.([ 

suggest that the tunneling effect 
is calculable in much the same way as for periodic potentials, and that in fact the 
result can be obtained by identification of appropriate parameters. We do not enter 
into an extensive calculation here along the lines of ref.||16||, and instead proceed 



by identification. One avoids ugly integrals by setting in the Schrodinger equation 
E = E^ +1 + A, where A is the perturbation theory correction of the eigenvalue. 
The original Schrodinger equation then becomes 








(11) 
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where 



Setting 



G 2 (0) 



2) 

'sin <b — a) 2 



(12) 



= H(0)exp j± J^G((p)d^ 



we obtain the WKB exponential 



exp I ± / G((j))d(f) \ = exp -j ± S (cos + a 



>/5 



(13) 



(note that without inclusion of the approximate form of E the integral would not 
have been so simple !). One can compute 5(0) etc. as in ref.[16]. The boundary 
conditions require the evaluation of the wave function above the chosen minimum of 
the potential (say that at 0o,i), implying for the barriers to the left and to the right 
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These expressions are seen to be (cf.ref. [|T^]) precisely the values of the action of the 
instantons travelling through the two differently sized barriers between (0o,i,0o,o) 
and (0o,i,0o,2) respectively as shown in Fig. |I] (one could say that these instantons 
describe clockwise and anticlockwise underbarrier rotations of the magnetic moment). 
In order to obtain this result it was essential to insert the oscillator approximated 
energy eq.fllOD into the Schrodinger equation. It is exponential factors like those of 
eq. ([i~3D with the boundary conditions of eq.flllj) which are typical of tunneling con- 
tributions. In the present case both of these contribute to the overall level splitting. 
Knowing these factors we can write down the level splitting by making the appro- 
priate replacements in the formula for the level splitting in the case of the Mathieu 
equation (cf. ref.|l6j ,also cited in [fl2|| ) and adding these with equal weights so that 
in the limit a->0we regain the level splitting of the case without the magnetic field. 
The factors multiplying the exponentials are characteristic of the central well (i.e. 

of eq.(^)) and in fact result from the matching of different branches of the wave 
function in domains of overlap as can be seen from ref . QT6J ; classically this factor de- 
scribes the number of bounces of the particle between the barriers before it escapes. 
Thus these factors are the same in both cases so that from the level splitting result 
of ref. [|l(J (also cited in |12|]) the level splitting of the nth excited state in the present 
case is obtained by the replacement 
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Figure 1: The periodic potential with asymmetric twin barriers and the instanton 
trajectories. 



and hence is found to be 
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In the limit a —>■ this reduces to the formula obtained in ref. |12| or to AEQ USt of 
formula (9a) of ref. . In particular we obtain for n = 



Ax 



16B 



g + 7 
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- (Vl— a 2 +a arcsin a) 



cosh(- 



(i±j) 



7T 



(16) 



which is to be compared with the corresponding path integral result of ref. || (there 
eq.(16)). In our result the origin of every factor is clearly understood as explained 
above - of course, with the assumption B > A and s > 1. The somewhat different 
factors in ref. M result from the complicated path integral calculation, taking into 
account of the field dependence of the mass and appropriate approximations. 

In Fig. we plot our expression A x for the values given in Fig. 2 of ref. ||. 
We make the amazing observation that the plots are practically identical in spite of 
the fact that the value of b chosen (i.e. b = 1) is really too small and still better 
agreement with the exact values will be obtained for larger values of b. 
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In FigsJ^ and |4] we display further plots of the level splitting demonstrating its 
increase with the magnetic field (as desired for better observability) and its variation 
with b for a fixed value of s. 

In Table |] we display some absolute values of the level splitting as calculated 
from our result and compare these with values given in ref . || . 

The very simple derivation of the nontrivial level splitting of eq.(|16D given here 
demonstrates the calculational superiority of the Schrodinger method in quantum 
mechanics. Of course, the path integral method has considerable pedagogical value 
in quantum mechanics as a model for more complicated applications in field theory. 
The derivation given here also shows clearly that the parameters s and have 
to be large for the asymptotic expansions of the solutions and eigenvalues to make 
sense. One can also see that the effect of the field dependence of the mass - here 
neglected - cannot be large (in agreement with our findings in ref.[^|). A further 
considerable advantage of the derivation given here is - apart from checking the 
path integral result - that it immediately supplies the splitting of higher degenerate 
oscillator levels which cannot be obtained with vacuum instanton methods. It is 
interesting to observe that the WKB exponential which corresponds to the classical 
or tree approximation is the same also for excited states. Thus if this is calculated for 
nonvacuum (or periodic) instantons the same expression must be obtained. We did 
not specify above whether s is integral or half-integral. The reason is that for half 
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Figure 3: The ratio of the level splitting divided by that for h = (i.e. magnetic 
field zero) plotted against a. One should note that a is effectively h divided by the 
large spin s. 




Figure 4: The level splitting at fixed value s = 10 for b = 5, 10, 15 and 20. In our 
considerations here b has to be larger than I. 
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h 


Ax 


Ai/AEo 


AE l nst /AE 


0.0 


2.4488 x 10" 4 


1.104 


1.029 


0.4 


3.4375 x 10" 4 


1.110 


1.065 


0.8 


7.0830 x 10" 4 


1.111 


1.123 


1.2 


1.5803 x lO" 3 


1.145 


1.174 



Table 1: Level splitting values calculated from Ai in units of A compared with the 
numerical values AE and the semiclassical results AE^ nst of ref. || for s = 10, b = 2 
and different values of h 



integer spins the classical twofold degeneracy remains due to Kramer's degeneracy 
||. For completeness we mention that the case B = has been dealt with in ref. 



the method there employed being one applicable to the original discrete spin system. 
The method discussed here obviously also shows how the level splitting in the case 
of a periodic potential with periodically recurring asymmetric twin barriers can be 
obtained. We remark that numerous properties of spin systems have been considered 
in QT9 1 . Finally our presentation should also appeal to experimentalists in view of 



the transparency of the calculational steps. 
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